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I.  INTRODUCTION 

A  height-of-burst  (HOB)  test  such  as  Pre-DIRECT  COURSE  creates  more  areas 
of  interest  than  simply  shock  phenomenology.  To  obtain  a  known  HOB  experiment 
with  a  high-explosive  charge  some  type  of  support  structure  is  required.  This 
support  structure  may  cause  hazards  from  fragments.  The  fragments  are  a 
result  of  the  large  intensity  blast  and  thermal  loading  of  the  support  struc¬ 
ture.  Given  a  certain  set  of  initial  conditions  a  fragment  could  be  capable 
of  traveling  tens  of  thousands  of  meters. 

Many  locations  in  the  test  area  require  that  personnel  be  present,  operat¬ 
ing  equipment  or  conducting  the  test  (i.e.,  instrumentation  vans,  administra¬ 
tion  vans,  laser  bunker,  weather  station,  and  observation  area).  The  possible 
presence  of  fragments  then  creates  safety  hazards  and  individual  experiment 
hazards.  Possibilities  will  exist  for  gages  or  objects  to  receive  damage  from 
fragment  impact. 

This  effort  will  attempt  to  predict  the  quantity,  velocity,  and  range  of 
fragments  present.  To  verify  the  predictions,  high  resolution  photographic 
coverage  will  be  utilized  during  the  event  from  both  ground  surface  and  aerial 
platforms.  The  final  phase  of  the  effort  will  consist  of  a  posttest  ground 
search  and  documentation  of  the  recovered  fragments. 
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II.  DESIGN  OF  THE  EXPERIMENT 

Design  of  Support  Structure--To  support  the  2.177  x  10**  kg  (24  tons)  of  explo¬ 
sives  at  a  height  of  17.37  m  (57  ft)  an  extensive  design  effort  was  initiated. 
Considering  all  of  the  design  factors,  (i.e.,  cost,  fidelity,  constructability, 
safety,  etc.)  a  single,  guyed,  battened  column  was  chosen  as  the  best  concept 
(Fig.  1).  This  column  extends  from  the  ground  surface  to  a  point  just  below 
the  charge  surface.  Attached  to  the  top  of  the  battened  column  is  a 
fiberglass  column.  The  fiberglass  column  passes  through  the  sphere  of  explo¬ 
sives  and  extends  1.29  m  (4.25  ft)  above  the  explosives.  On  top  of  the 
fiberglass  column  is  an  upper  steel  column  2.44  m  (8  ft)  long. 

The  explosives  are  held,  with  the  center  at  17.37  m,  in  a  spherical 
geometry  by  a  fiberglass  container.  This  container  is  supported  in  place  by  a 
polyester  webbed  net  which,  in  turn,  is  connected  to  a  series  of  wire  ropes 
which  transfer  the  load  to  the  top  of  the  upper  steel  column  (Fig.  2). 

Fragments--Al  1  of  the  members  making  up  the  support  structure  will  create  frag¬ 
ments  of  some  type,  but  not  all  of  them  will  create  safety  hazards.  The  upper 
portion  of  the  lower  battened  column  will  receive  tremendous  loadings  from  the 
blast  and  thermal  environments.  As  was  seen  at  the  preliminary  test,  FRED  II, 
the  upper  one-third  of  the  battened  steel  column  was  stripped  and  scattered 
(Ref.  1).  However,  these  fragments  are  scattered  at  small  ranges  from  ground 
zero.  This  part  of  the  support  structure  will  not  cause  safety  hazards  since 
the  fragments  will  be  directed  downward  due  to  their  initial  location,  under¬ 
neath  the  charge.  The  maximum  ranges  from  these  fragments  will  probably  be 
gained  from  bounces  (after  the  fragment  has  initially  struck  the  concrete 
radial s).  The  major  hazards  these  fragments  would  produce  would  be  to  various 
experiments  in  the  region  from  ground  zero  to  approximately  two  HOBs.  These 
fragments  when  initially  impacting  the  ground  could  have  velocities  up  to 
hundreds  of  meters  per  second  and,  thus,  create  substantial  damage. 

The  fiberglass  column  and  container  will  produce  fragments,  but  will  not 
create  substantial  safety  or  experiment  hazards.  These  fragments  may  obtain 
velocities  up  to  thousands  of  meters  per  second,  but  due  to  their  low  mass 
will  not  carry  much  momentum  (relative  to  the  steel  members).  The  major 
hazards  would  be  to  transmission  experiments  or  direct  measurement  devices. 
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Figure  1.  Pre-DIRECT  COURSE  concept. 
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Figure  2.  Upper  portion  of  support  structure. 
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The  remaining  support  structure  member  (the  upper  steel  column)  has  the 
greatest  potential  for  creating  safety  hazards.  Included  in  the  fragments 
from  this  column  are  the  hooks,  rings,  and  bolts  which  connect  the  polyester 
net,  wire  ropes,  and  guy  cables.  Since  all  of  these  objects  are  located  above 
the  charge  surface,  their  potential  for  long  flights  is  the  greatest  of  any 
member  of  the  support  structure. 

Upper  Steel  Column--The  upper  steel  column  is  essentially  a  standard  40.64  cm 
(16  in)  outer  diameter  steel  pipe  2.439  m  (8  ft)  long  with  0.95  cm  (3/8  in) 
walls.  Welded  to  the  pipe  are  various  flanges  and  brackets  for  connecting  guy 
and  wire  rope  cables  (Fig.  3). 

The  pipe  is  ASTM-53A  steel  having  a  yield  strength  of  2.482  x  10®  Pa 
(36  ksi),  which  represents  approximately  60  percent  of  the  ultimate  strength. 
Flanges  and  plates  are  also  ASTM-53A  steel,  and  are  attached  with  welds  with  a 
strength  of  4.826  x  10®  Pa  (70  ksi)  making,  in  essence,  a  stronger  material  at 
connection  points  than  the  steel  material  itself. 

Method  of  Fragmentation  —  Due  to  the  design  of  the  center  fiberglass  column,  it 
is  assumed  that  the  interior  of  the  column  will  act  as  a  shock  tube  confining 
a  portion  of  the  upward  flow  to  two  dimensions  or  less.  The  flow  outside  of 
the  column  will  be  fully  three-dimensional  and  diverge  spherically.  This  dif¬ 
ference  would  allow  the  flow  inside  the  column  to  decelerate  slower  than  the 
exterior  flow,  if  friction  and  other  small  losses  are  neglected.  Assume  that 
the  shocked  flow  has  propagated  for  some  amount  of  time  and  reached  the  upper 
steel  column.  The  shock  which  has  been  confined  by  the  shock  tube  effect 
would  be  loading  the  interior  of  the  pipe  prior  to  the  loading  of  the  pipe's 
exterior.  This  pressure  differential  could  result  in  a  loading  gradient  which 
could  stress  the  pipe  section  to  a  point  of  fracture.  Fortunately,  the 
majority  of  the  pipe  surface  is  parallel  to  the  direction  of  the  flow  and 
receives  most  of  the  loading  from  static  pressures  and  little  from  dynamic 
pressures.  The  flanges  and  plates,  however,  would  receive  a  larger  loading 
due  to  the  stagnation  of  the  flow. 

Many  problems  are  faced  when  trying  to  predict  the  fragmentation.  Also, 
the  environment  at  these  close  ranges  is  nearly  undefinable.  There  is  no 
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documented  evidence  of  recorded  pressures  or  impulses  at  ranges  of  1.37  m  to 
3.809  m  (4.5  ft  to  12.5  ft)  from  the  surface  of  a  large,  high  explosive 
charge.  Photographic  data  may  be  obtained  relating  to  the  temperatures  and 
velocities,  but  these  data,  at  the  ranges  of  interest,  cannot  be  used  to 
deduce  pressure  or  density  information.  Normal  shock  relations  are  not  valid 
since  the  shock  and  fireball  are  still  attached.  The  only  means  to  obtain 
information  is  from  theoretical  or  hydrocode  calculations.  But,  once  again, 
at  these  regions,  these  numerical  data  have  not  been  validated. 

There  is  no  documented  evidence  of  any  method  to  predict  random  fragmen¬ 
tation  of  any  material.  Hydrocodes,  such  as  Elastic-Plastic  HULL  (Ref.  2), 
are  available  which  can  model  the  steel  members  and  the  explosive  material. 
However,  the  most  these  calculations  will  reveal  is  the  maximum  stress  and 
strains,  and  the  energy  transfer  into  the  materials,  and  nothing  about 
fragment  shapes  or  sizes.  This  information,  about  the  strain  and  stresses, 
can  be  useful  to  determine  if  the  member  will  be  loaded  beyond  its  ultimate 
strength. 

Conventional  weapon  research  does  reveal  some  information  concerning  frag¬ 
ment  trajectories  if  the  explosive  weight,  casing  parameters,  and  materials 
are  known.  Once  again,  these  methods  reveal  little  information  about  fragment 
sizes  or  shapes.  Another  problem  with  using  conventional  weaponry  research  is 
the  absence  of  any  data  with  low  explosive  characteristics  (i.e.,  detonation 
pressures  and  velocities)  such  as  Ammonium  Nitrate  and  Fuel  Oil  (ANFO). 

Past  large-scale  testing  and  results  from  any  conventional  measurement 
schemes  do  not  reveal  any  information  on  break-up  of  close-in  structures.  The 
structures  in  these  close-in  regions  are  engulfed  by  the  fireball  before  any 
noticeable  failure  or  movement  occurs.  There  are  two  methods  by  which  some 
data  could  be  gained  which  are  infrared  photography  and  Doppler  radar 
transmission.  However,  the  resolution  of  the  infrared  photography  is  too  poor 
to  give  adequate  break-up  data  and  the  radar  transmission  is  hard  to  field  and 
very  expensive. 

Given  these  obstacles  an  attempt  was  made  to  use  bits  and  pieces  of  the 
various  methods  to  predict  the  fragment  parameters. 
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Calculation  Methods — As  previously  mentioned  there  is  no  evidence  of  one 
method  to  cover  the  entire  scope  of  this  effort.  Analytical  and  empirical 
methods  are  available  for  predicting  high-explosive  environments.  Methods  are 
available  for  predicting  environment  and  target  interactions.  Other  methods 
also  exist  for  predicting  fragment  trajectories.  However,  none  of  these 
methods  have  been  coupled  together  nor  can  they  predict  a  material's 
dissassembly  or  fragmentation. 

The  HULL  hydrodynamics  code  has  gained  international  recognition  for  its 
fairly  accurate  results  in  predicting  high-explosive  environments  (Ref.  3). 
S-CUBED  conducted  HULL  calculations  of  the  Pre-DIRECT  COURSE  event.  Data  from 
the  calculation  at  1.5,  2,  and  3  m  above  the  charge  surface  are  presented  in 
Figures  4a  through  4o  and  peak  values  are  summarized  in  Table  1  (Ref.  4). 

These  data  reflect  semifree-field  data.  Specifically,  the  plots  indicate  the 
environments  in  a  modeled  pipe  above  the  charge  surface.  The  pipe  is  non¬ 
reactive  with  the  environment,  with  the  exception  of  being  able  to  reflect  the 
flow  off  its  walls.  These  data  were  received  too  late  to  be  included  in  the 
predictions  for  this  effort,  but  are  presented  solely  for  completeness. 

One  hydrodynamic  method  which  was  utilized  as  a  source  of  initial  condi¬ 
tions  was  the  Elastic-Plastic  HULL  Code.  This  calculation,  performed  at  the 
Air  Force  Weapons  Laboratory  (AFWL),  modeled  the  explosive  charge  and  the  upper 
support  column,  complete  with  flanges.  The  fiberglass  around  the  charge  and 
the  center  column  was  not  modeled.  This  calculation,  problem  No.  40040.0010, 
consisted  of  a  cylindrical  coordinate  system  with  close-in  cell  dimensions  of 
5  cm  radially  by  5  mm  vertically.  The  explosive  modeled  was  ANFO  and  the 
steel  was  ASTM-53A. 

Figure  5  shows  a  density  plot  at  500  pS.  The  upper  support  column  and  the 
shock/fireball  location  can  be  identified  (the  mesh  is  originated  at  the 
center  of  the  charge).  At  this  point  in  calculation  time  the  detonation  wave 
has  propagated  through  the  charge  and  the  shock/fireball  has  grown  by  80  cm. 
Figures  6  through  10  show  additional  density  and  velocity  vector  plots  from 
the  calculation.  One  of  the  many  drawbacks  to  hydrodynamic  calculations  is 
their  enormous  cost;  for  this  reason  this  calculation  was  ended  near  1  ms.  At 
this  point  the  calculation  indicated  that  the  center  of  mass  of  the  steel  had 
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Figure  4.  Calculation  results. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  4.  Continued. 
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Figure  5.  E-P  HULL  density  results  at  500  ys. 
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Figure  6.  E-P  HULL  vector  velocity  results  at  500  ys. 
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Figure  7.  E-P  HULL  density  results  at  736  ys. 
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Figure  8. 


E-P  HULL  vector  velocity  results  at  736  ys. 
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Figure  9. 


E-P  HULL  density  results  at  948  ys, 
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E-P  HULL  vector  velocity  results  at  948  ys. 
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a  radial  velocity  of  60.1  m/s  (197.2  ft/s),  an  axial  velocity  of  61.9  m/s 
(203.1  ft/s),  (or  a  speed  of  86.3  m/s  (283.2  ft/s))  directed  upward  from  the 
vertical  axis  at  an  angle  of  45.85  deg.  Since  the  calculation  can  relate  the 
transfer  of  energy  from  the  environment  to  the  target  material,  and  does 
account  for  the  strength  properties  of  the  target,  some  deformation  is 
observed.  However,  the  deformation  at  the  end  of  this  calculation  indicates 
that  the  pipe  is  stretching  radially,  but  has  not  broken.  Information  can  not 
be  gained  about  individual  fragment  geometries  or  masses.  Another  important 
concern  is  the  time  duration  of  the  calculation.  Due  to  the  lack  of  computa¬ 
tional  time,  the  shock/fireball  has  not  totally  engulfed  the  upper  support 
column.  This,  in  turn,  means  that  the  total  impulse  associated  with  the  event 
has  not  been  imparted  to  the  upper  structure;  therefore,  the  mass,  velocity 
and  angle  data,  previously  mentioned,  are  preliminary.  The  data  are  for  the 
material's  parameters  after  the  start  of  the  loading  process.  In  actuality 
these  data  represent  a  point  along  a  dynamic  curve  which  has  not  yet  reached 
steady-state  conditions  (i.e.,  the  upper  steel  column  is  still  accelerating 
and  deforming). 

Conventional  weaponry  research  has  revealed  information  about  initial  con¬ 
ditions  of  bomb  fragments.  These  conditions  are  deduced  from  measurements 
when  certain  constraints  are  placed  on  the  test  bomb.  Usually  the  bomb  is 
made  up  of  highly  controlled  explosives  and  the  geometry  is  cylindrical  or 
spherical.  The  casing  of  the  explosive  is  controlled  by  defining  a  desired  . 
geometry  by  scoring.  Through  high-speed  photography  and  impact  screens, 
fragment  velocities  and  flight  characteristics  can  be  recorded.  Knowing  the 
aerodynamic  characteristics  of  the  known  fragment,  the  initial  velocity  can  be 
backed  out.  Knowing  the  initial  velocity,  the  energy  imparted  to  the  fragment 
from  the  known  explosive  can  be  deduced.  The  energy  transferred  to  the 
fragment  is  known  as  the  Gurney  energy  constant,  or  the  energy  per  unit  mass 
of  the  explosive  that  can  be  converted  into  mechanical  work  (Ref.  5).  Using 
the  Gurney  energy  constant,  bomb  fragment  distributions  can  be  calculated. 
These  distributions  are  most  easily  calculated  using  the  work  of  Bishop  (Ref. 
6).  The  biggest  problem  is  associating  this  type  of  fragment  analysis  with  a 
low  energy  explosive  such  as  ANFO.  Another  problem,  when  extending  conven¬ 
tional  weaponry  work  to  this  specific  event,  is  the  location  of  the  fragment 
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material.  In  conventional  weapons  the  metal  material  surrounds  the  explosive, 
however,  in  this  effort  the  metal  material  is  located  almost  one  charge  radius 
from  the  explosive  surface. 

In  an  attempt  to  relate  known  Gurney  energy  constants  to  explosives,  it 
was  found  that,  as  suspected,  the  Chapman-Jougnet  (C-J)  pressure  is  directly 
related  to  the  energy  constants.  Figure  11  shows  a  plot  of  15  Gurney  energy 
constants  (Ref.  7)  (/^)  versus  C-J  detonation  pressures  (Ref.  8).  It  was 
assumed  that  zero  energy  would  produce  zero  C-J  pressure  and,  therefore,  a 
zero  value  of  Gurney  energy  constant.  Thus,  a  polynominal  of  the  second  order 
provided  a  close  fit  to  these  data.  The  form  of  the  polynominal  is  as 
fol  1  ows : 

(/2E  X  10^)  =  2.20214  -  45.9372  x  -5.43194  x^ 

2 

where  x  =  C-J  pressure  x  10  kbar. 

Using  the  value  of  61  kbar*  as  the  C-J  detonation  pressure  for  ANFO  a 
Gurney  energy  constant  of  854  m/s  (2.802  ft/s)  can  be  calculated.  This  value, 
while  not  validated,  will  be  assumed  to  be  adequate  for  the  ANFO. 

To  summarize,  S-CUBED  calculated  that  the  pipe  will  receive  a  free-field 
pressure  on  the  order  of  100  MPa.  The  Elastic-Plastic  HULL  calculation 
revealed  that,  at  near  1  ms,  the  speed  of  the  center  of  mass  of  the  steel  was 
86.3  m/s.  Finally,  from  conventional  weaponry  work,  the  Gurney  energy 
constant  is  approximated  as  854  m/s. 

Fragment  Trajectories--Many  methods  are  available  for  trajectory  analysis, 
most  of  which  stem  from  aerodynamics.  The  simplest  form  of  trajectory  analy¬ 
sis  is  presented  in  Appendix  A.  This  derivation  of  trajectory  assumes  that 
the  drag  force  is  small  relative  to  the  gravity  force.  This  assumption  is 
adequate  for  ballistic,  streamlined  objects  but  may  not  represent  the  frag¬ 
ments  of  concern  in  this  effort.  The  range  equation,  derived  from  this  analy¬ 
sis  is  as  follows: 

Vi^ 

Rc  =  — _  si n  $,• 

g 


*Private  communication  with  C.E.  Needham,  S-CUBED,  Albuquerque  NM,  February 
1983. 
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Figure  11.  Gurney  energy  constants  versus  C-J  pressures. 
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where 

Rg  =  maximum  range 

V^-  =  initial  velocity 

g  =  acceleration  due  to  gravity 

=  initial  angle  of  projection 

The  major  downfall  of  this  analysis  is  the  assumption  that  the  force  due  to 
drag  is  negligible  relative  to  the  force  due  to  gravity.  The  flight  of  any 
fragment  in  Pre-DIRECT  COURSE  is  relatively  close  to  the  surface  of  the  earth 
so  the  density  of  air  also  becomes  important.  Since  the  drag  coefficient  is 
large  (one  or  greater),  it  cannot  be  assumed  that  the  fragment's  drag  coef¬ 
ficient  is  low  (i.e.,  the  fragment  does  not  happen  to  be  streamlined). 

Detailed  conventional  weapon  analysis,  such  as  derived  in  Reference  6 
assumes  that  drag  is  a  major  factor.  This  analysis  accounts  for  fragment  size 
and  shape  as  well  as  the  weight  of  the  explosive.  The  fragment  size  (or 
thickness)  is  resolved  to  be  the  most  important  factor  in  determining  the 
maximum  fragment  range.  It  is  also  stated  in  Reference  6  that  the  greatest 
possible  fragment  velocity  occurs  when  the  case  weight  is  negligible  compared 
to  the  explosive  weight.  In  Pre-DIRECT  COURSE  the  upper  support  column  has  a 
mass  of  approximately  453.6  kg.  Thus,  if  it  is  assumed  that  this  mass  is 
evenly  distributed  around  the  exterior  of  the  explosive,  the  case  to  the 
explosive  weight  ratio  is  1/46  (a  relatively  small  number  compared  to  conven¬ 
tional  weapons). 

Bishop  (Ref.  6)  utilizes  the  aerodynamic  drag  as  established  by  Hoerner 
for  regular  shapes  in  subsonic  and  supersonic  velocities.  Two  distinct 
classes  of  shapes  are  considered.  The  first  class  is  composed  of  regularly 
shaped  projectiles  with  sharp  pointed  conical  or  wedge  shaped  noses.  However, 
due  to  probability  analysis  it  is  unlikely  that  this  would  occur.  Assume, 
though,  that  one  fragment  did  possess  these  properties.  Generally  this 
fragment  would  not  have  its  center  of  mass  located  to  insure  stable  motion 
through  the  air  with  the  sharp  point  or  edge  headed  forward  at  all  times. 
Bishop's  second  class  of  shapes  are  those  that  do  not  have  the  special  charac¬ 
teristics  of  the  first  class.  These  fragments  are  characterized  by  a  constant 
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coefficient  of  drag  over  a  wide  range  of  supersonic  velocities.  The  coef¬ 
ficient  of  drag  is  always  greater  than  one  except  for  regularly  shaped  rounded 
bodies  (i.e.,  spheres  and  cylinders)  for  all  sub-  and  supersonic  velocities. 

The  coefficient  remains  as  such  at  trans-  and  supersonic  velocities  due  to  the 
occurrence  of  a  detached  shock  wave  some  distance  in  front  of  the  projectile. 

The  shock  wave  is  detached  for  all  fragments  except  those  in  the  first  class. 

In  the  first  class  the  shock  wave  can  be  attached  to  the  nose  of  the  fragment 
only  when  the  nose  angle  is  less  than  e,  where  e  =  1/M  and  M  =  the  Mach  number. 
All  other  fragments  are  in  the  second  class  and  have  a  detached  shock  wave. 

Finally,  consider  the  orientation  of  the  fragment  in  flight.  The  general 
condition  for  approximately  stable  flight  of  a  nonrotating  fragment  requires 
that  the  center  of  mass  be  well  ahead  of  the  line  of  action  of  the  resultant 
aerodynamic  drag  force  for  at  least  a  small  range  of  fragment  orientations  in 
any  direction  from  a  zero-lift  attitude.  Otherwise  the  particular  zero-lift 
attitude  assumed  is  unstable,  and  the  fragment  tends  to  oscillate,  or  even 
tumble.  A  homogeneous  nonrotating  fragment  having  a  uniform  density  and  an 
unstable  zero-lift  orientation  is  usually  reached  when  the  fragment  presents 
its  minimum  frontal  projected  area.  This  applies  especially  to  a  thin  flat 
sheet  of  material  (but  not  rotating). 

A  thin  flat  sheet  is  generally  unstable  when  moving  edge-on  through  the 
air,  since  the  aerodynamic  lift  force  (at  any  small  angle  of  attack)  acts  with 
a  torque  couple  about  the  center  of  mass--which  tends  to  rotate  the  sheet  into 
a  position  with  a  larger  angle  of  attack.  The  larger  angle  of  attack  results 
in  an  increased  torque,  which  tends  to  increase  the  angle  of  attack  even 
further,  etc.  Thus,  the  leading-edge  type  of  flight  is  an  unstable  zero-lift 
orientation. 

The  perpendicular  type  of  zero-lift  orientation  is  at  least  semistable; 
when  a  flat  shape  moves  so  that  the  maximum  frontal  area  is  presented,  there 
is  a  range  of  the  angle  of  inclination  within  which  the  resulting  torques  are 
of  a  restoring  variety  instead  of  an  overturning  variety.  The  restoring 
torque  (approximately  cosine)  is  small,  which  means  that  the  flat  shape  may 
oscillate  through  a  considerable  angle  before  any  appreciable  restoring  torque 
comes  into  play.  Thus,  the  slightest  air  flow  disturbance  results  in  direc¬ 
tional  oscillations  which  do  not  damp  out  readily.  When  the  fragment  turns 
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completely  over,  instead  of  merely  oscillating  back  and  forth,  a  more  compli¬ 
cated  analysis  is  required.  This  analysis  is  presented  in  Reference  6  work 
and  will  not  be  repeated  here. 

It  is  stated  in  Reference  6:  "One  can  be  sure  that  the  maximum  dimensions 
of  any  fragment  will  be  less  than  the  radius  of  the  container  surrounding  the 
explosive  charge."  However,  when  an  explosive  charge  is  not  closely  contained 
or  when  asymmetries  exist,  the  fragment  may  be  large,  as  in  Pre-DIRECT  COURSE. 
However,  Bishop  does  not  account  for  fragments  considered  in  the  specific  case 
of  the  Pre-DIRECT  COURSE.  He  does  state  that  the  edge  and  face  area  results 
in  an  aerodynamic  factor  that  is  "a  statistical  quantity  that  can  only  be 
found  by  observing  a  large  number  of  fragments  from  many  experiments". 

Even  though  Gurney's  and  Bishop's  analyses  do  not  consider  the  specific 
case  of  Pre-DIRECT  COURSE,  ball  park  factors  may  be  gained  from  this  type  of 
analysis.  In  the  conventional  exercises  the  majority  of  the  explosive  energy 
is  directly  imparted  to  the  encasing  metal.  In  Pre-DIRECT  COURSE,  the  metal  is 
almost  one  charge  radius  away  from  the  explosive.  When  considering  spherical 
divergence,  it  can  be  safely  assumed,  that  less  energy  is  imparted  to  the 
metal  than  if  it  surrounded  the  explosives. 

A  summary  of  the  equations  derived  by  Bishop  are  presented  in  Appendix  B. 
This  analysis  will  be  used  in  predicting  fragment  trajectories  for  Pre-DIRECT 
COURSE.  The  values  from  this  analysis  are  assumed  to  be  liberal  considering 
the  aforementioned  caveats. 

Fragment  Trajectory  Calculations--Table  2  presents  the  presumed  fragment  ini¬ 
tial  conditions  from  the  calculations  and  extrapolated  Gurney  analysis.  These 
values  represent  the  Gurney  energy — density  constants  only  and  do  not  factor 
in  the  geometrical  shape  of  the  bomb.  From  Bishop's  analysis,  one  realizes 
that  as  the  ratio  of  explosive  weight  to  case  weight  increases,  the  initial 
velocity  of  the  fragments  (as  predicted  by  Gurney)  approaches  the  value  of 
/^.  Since  Pre-DIRECT  COURSE  had  a  low  weight  ratio  it  can  be  assumed  that 
the  predicated  and  calculated  values  of  /2E  are  the  initial  velocities. 

Utilizing  these  initial  conditions  the  various  ranges  were  calculated 
using  the  methods  derived  and  explained  in  Appendices  A  and  B.  The  results 
from  these  exercises  are  presented  in  Appendix  C.  To  estimate  fragment  sizes. 


35 


AFWL-TR-83-52 


TABLE  2.  FRAGMENT  INITIAL  CONDITIONS 


Method 


Initial  Velocity 


Gurney-Energy 
Density  Constant 


S-CUBED  Calculations 


UMDEF 


NA 


Elastic-Plastic  HULL  Calculations  86.3  n/s  at  1  ms 


MA 


C-J  Pressure  vs  Analysis 


854  m/s 


2,802  ft/s 


in  the  Bishop  analysis  case,  random  dimensions  were  chosen  to  try  to  bound  the 
expected  sizes.  Figure  12  presents  the  fragment  size  versus  calculated  range 
for  a  given  initial  velocity.  These  values  must  be  used  with  caution  since, 
as  previously  stated,  the  simplified  analysis  is  a  known  overestimation  and 
Bishop's  analysis  is  indirectly  related  to  Pre-DIRECT  COURSE. 

Experimental  Procedures-- In  an  attempt  to  verify  the  theoretical  calculations 
associated  with  this  effort,  high  resolution  photographic  measurements  are 
required.  This  coverage  requires  both  ground  and  aerial  platform-based 
coverage.  The  photography  will  consist  of  both  70  mm  and  35  mm  ground-based 
and  35  mm  aerial  photography.  The  cameras  will  have  overlapping  formats  and 
record  at  speeds  of  30  and  360  fr/s.  Three  ground-based  camera  stations  will 
be  used  in  an  attempt  to  correlate  one  camera's  images  with  another.  The  data 
from  these  cameras,  if  successful,  will  give  tighter  and  more  precise  data 
than  from  any  one  camera.  One  problem  with  noncorrel ated  ground-based  photo¬ 
graphy  is  the  image  plane  of  the  camera.  A  camera  will  record  anything  in  the 
field  of  view.  This  field  includes  any  objects  moving  toward  or  180  deg  away 
from  or  anywhere  in  between.  To  precisely  calculate  fragment  trajectories, 
the  fragment  must  be  on  a  plane  perpendicular  to  the  camera's  head  on  viewing 
angle.  If  the  fragment  is  traveling  on  any  other  course  the  cosine  of  the 
angle,  relative  to  the  perpendicular  must  be  considered. 

To  solve  this  problem,  two  options  become  available.  These  are:  (1)  an 
increased  quantity  of  ground-based  cameras,  or  (2)  aerial  cameras  which  pro¬ 
vide  another  plane  of  the  three-dimensional  geometry.  Therefore,  driven  by  a 
cost  analysis  basis  it  would  be  advantageous  to  utilize  aerial  photography. 
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Figure  12.  Fragment  size  versus  ground  range. 
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In  addition  to  the  dynamic  photography,  two  additional  validation  tech¬ 
niques  are  to  be  utilized.  The  first  is  aerial  photography.  This  photography 
would  be  low  level,  high  resolution  ground  mapping  of  the  surface.  Pretest 
aerial  photography  will  be  taken  of  the  test  bed  area  out  to  609.5  m  radius 

from  ground  zero.  Posttest  photography  will  repeat  this  documentation. 

Any  differences  in  pre-  and  posttest  documentation  will  be  recorded  and  be 
investigated.  The  second  validation  technique  is  a  rigorous  ground  search. 

This  ground  search  will  consist  of  an  extensive  walk  search  for  fragments  from 
ground  zero  out  to  ranges  of  609.5  m.  Once  the  fragments  are  tagged,  another 
team  will  survey  and  record  the  specific  fragment  data. 

The  coupling  of  these  measurement  techniques  is  vital  to  verifying  the 
predictive  methodologies.  If  any  portion  of  the  technique  is  omitted  much 
more  analysis  is  required  and  will  not  yield  as  accurate  results. 

Structural  Identification--To  identify  the  fragments  during  the  posttest  recov¬ 
ery,  several  steps  were  taken  to  uniquely  mark  various  pieces.  The  top  26 
channels  of  the  battened  column  were  indented,  by  metal  stamping,  with  a 
letter  of  the  alphabet.  The  letter  A  was  on  the  top  rung,  progressing  down, 
rung  by  rung  to  the  letter  Z.  In  addition  to  the  letter,  each  side  was  stamped 
with  a  number  (1  through  4).  The  sequence  of  number  and  letter  was  stamped 
several  times  on  each  channel,  in  anticipation  that  the  channel  would  break. 
This  method  of  identification  covered  each  side  of  the  battened  column  for 
more  than  50  percent  of  the  length  of  the  entire  column. 

The  upper  steel  column  was  also  marked.  The  pipe  was  also  metal  stamped 
using  the  aforementioned  procedure.  The  letter  A  started  at  the  bottom  of  the 
pipe  (next  to  the  attachment  to  the  fiberglass  column)  and  extended  by  10.16  cm 
increments,  to  the  letter  Y.  The  numbers  1  through  4  were  again  used  with  the 
letters  to  form  four  separate  vertical  lines  (90  deg  apart)  of  identification 
over  the  entire  length  of  the  pipe.  In  addition  to  the  markings,  the  upper 
support  column  and  attachment  were  painted  fire  engine  red. 

The  numbered  vertical  lines  on  the  upper  support  column  and  the  numbered 
sides  of  the  battened  column  were  lined  up.  The  number  1  line  was  oriented 
toward  33  deg  (test  bed  coordinates),  number  2  toward  123  deg,  number  3  toward 
213  deg,  and  number  4  toward  303  deg. 


38 


AFWL-TR-83-52 


Jet  Retardation--Sufficient  concern  was  generated  over  the  shock-tube  effect 
that  a  jet  retardation  scheme  was  designed.  The  goal  was  to  create  a  baffle 
inside  the  upper  portion  of  the  fiberglass  column  and  the  upper  support  column 
which  would  control  a  jet  and  allow  the  exterior  of  the  pipe  to  be  loaded 
first.  If  the  jet  could  be  controlled,  the  upper  pipe  would,  theoretically, 
be  put  into  compression  and  possibly  not  fragment  to  a  large  degree. 

Figure  13  shows  the  scheme  designed  for  Pre-DIRECT  COURSE.  The  system  is 
composed  of  polystyrene,  plywood  disks,  and  sand-filled  bags.  The  plywood 
disks  are  located  just  above  the  holes  in  the  fiberglass  column.  The  poly¬ 
styrene  disks  are  used  as  spacers.  The  attempt  is  to  retard  the  jet  enough  to 
fracture  the  fiberglass  column  near  the  area  of  the  holes.  In  this  manner. 


Figure  13.  Jet  retardation  system. 
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the  jet  may  be  damped  out  prior  to  reaching  the  steel  member.  The  purpose  of 
the  sand  bags  is  to  essentially  provide  more  mass  than  ambient  air  and 
possibly  delay  any  pressure  gradient  between  interior  and  exterior  loading  on 
the  pipe.  This  scheme,  although  primitive,  was  the  most  acceptable  due  to 
time  and  cost  limitations. 


Trajectory  Control --The  fragment  prediction  analysis  indicated  that  some  of 
the  upper  support  column  would  not  be  fragmented.  The  possibility  of 
controlling  the  trajectory  of  the  remaining  column  was  considered.  The  New 
Mexico  Engineering  Research  Institute  (NMERI)  designed  a  simple,  low  cost 
guidance  system.  Briefly,  cables  are  attached  at  one  end  to  the  upper  support 
column  and  the  other  end  is  attached  to  a  large  concrete  block.  The  block  is 
on  rollers  which  are  on  a  small  concrete  slab. 


As  the  upper 
concrete  block, 
jerk  the  cables, 
initially  direct 


support  column  moves,  the  cables  become  taut  and 
In  a  short  time  the  concrete  block  will  roll  off 
This  sudden  stop  will  pull  on  the  upper  support 
it  toward  the  concrete  block. 


pull  the 
the  slab  and 
column  and 
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III.  RESULTS 

Posttest  Overview — The  test  site  was  recovered  immediately  after  the  detona¬ 
tion.  From  a  distance,  a  large  portion  of  the  lower  battened  column  appeared 
to  be  standing,  but  was  bent  and  twisted.  Upon  closer  examination  it  was 
noted  that  the  top  portion  of  the  column  extended  to  8.2  m  (26.9  ft)  above  the 
ground  surface.  An  even  closer  examination  revealed  that  only  approximately 
1.5  m  (4.92  ft)  of  the  column  was  left  standing  (Fig.  14).  This  photograph 
shows  that  the  column  was  sheared  off  at  the  four  uprights  approximately  1.5  m 
above  the  ground  surface.  The  upper  8.2  m  of  the  column  was  driven  down  and 
wedged  into  the  bottom  1.5  m  (Fig.  15).  Thus,  the  majority  of  the  battened 
column  just  happened  to  be  standing. 

Figure  16  shows  the  top  remaining  portion  of  the  battened  column.  Note 
that  the  channels  were  stripped  off  of  the  angle  iron  legs.  The  letter  guide 
system  indicated  that  the  channels  lettered  L  remained.  This  indicates  that 
eleven  levels  of  channel  were  stripped  off  the  column.  Also  note  that  some  of 
the  angle  iron  legs  are  missing.  The  points  of  failure  of  the  channels  and 
angle  iron  appeared  to  be  at  locations  of  welds.  .  The  welds  held  but  the 
material  fractured  at  the  high  local  stress  locations. 

Figure  17  shows  a  close-up  of  the  middle  portion  of  the  battened  column. 

A  distinctive  bend  is  noted.  Correlation  back  to  pretest  conditions  indicates 
that  this  is  the  location  of  the  polystyrene  block.  Apparently  as  the 
shock/fireball  was  traveling  down  the  column  the  polystyrene  acted  as  a  tem¬ 
porary  blockage  stagnating  the  flow  and  causing  a  high  reflected  pressure 
which  deformed  the  column. 

The  typical  fragments  of  concern  in  this  study  are  shown  in  Figures  18  and 
19.  Figure  18  shows  a  piece  of  channel  of  approximately  105  cm^  (16.28  in^). 
The  final  location  of  this  fragment  was  97.5  m  (319.9  ft)  from  ground  zero. 
Figure  19  shows  a  larger  fragment  (the  dark  object  in  the  lower  left  corner). 
This  fragment  was  approximately  310  cm  (48.05  in  ).  It  is  readily  apparent 
that  some  damage  could  have  easily  occurred  to  the  gage  mounts  in  the 
background . 
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Fiqure  14.  Remaining  structure  overview. 


Figure  15.  Close-up  of  bottom  of  structure 
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Figure  16.  Close-up  of  top  of  structure. 


Figure  17.  Close-up  of  middle  of  structure. 
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Figure  18.  Typical  small  fragment. 


Figure  19.  Typical  large  fragments. 
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Figures  20  and  21  show  the  damage  inflicted  upon  the  hardened  concrete 
radial s.  The  concrete  was  pulverized  from  the  direct  impacts  and  scattered 
over  the  area.  Figure  22  displays  the  metal  which  was  removed  from  the  hole 
shown  in  the  center  of  the  photograph.  Apparently  the  lower  column  acted  as  a 
funnel  for  some  period  of  time  and  allowed  the  fragments  to  drive  themselves 
into  the  concrete.  The  final  fragment  was  found  80  cm  (2.625  ft)  below  the 
concrete  surface.  The  total  mass  of  the  metal  found  in  the  hole  was  36.3  kg. 

The  major  portion  of  the  upper  steel  column  is  shown  in  Figure  23.  Approx¬ 
imately  60  percent  of  the  original  column  is  shown.  This  portion  was  found  on 
the  323  deg  azimuth,  22  m  (72.2  ft)  from  ground  zero.  One  can  easily  detect 
that  the  column  was  peeled  back  until  shear  occurred. 

Photographic  Analysis — The  experimental  photographic  information  provided  ade¬ 
quate  data.  The  ground-based  photography  was  analyzed  by  the  Physical 
Sciences  Laboratory  (PSL).  The  PSL  found  that  the  35  mm  photography  covered 
too  wide  of  a  range  and  could  not  provide  sufficient  resolution  to  track 
fragments  of  interests.  The  70  mm  photography,  although  limited  in  field  of 
view,  provided  enough  resolution  to  track  a  fragment  of  85  cm  (13.2  in  ). 

From  the  three  70  mm  cameras,  92  fragments  were  tracked.  The  results  from 
this  analysis  are  presented  in  Appendix  D.  These  data  reflect  calculated 
fragment  position-time  data,  as  well  as  predicted  velocities  and  suspected 
trajectory  angles.  Table  3  presents  a  summary  of  the  ground-based  photographic 
data.  Listed  in  this  table  is  the  apparent  area  of  the  fragment,  the  flight 
profile,  and  the  average  velocity.  The  minimum  and  maximum  values  from  this 
table  are:  apparent  area:  min — 0.0085  m^  (0.0915  ft^),  max— 0.4949  m^ 

(5.33  ft^),  average  velocity:  min— 32  m/s  (105  ft/s),  max — 629  m/s  (2064  ft/s). 
Also,  the  majority  of  the  observed  fragments  appeared  to  have  been  tumbling 
during  their  trajectories. 

Figures  24  through  26  depict  the  individual  fragment  trajectories  as 
viewed  from  the  three  camera  stations.  Note  that  a  large  circular  area  around 
the  depicted  charge  is  void  of  any  fragments.  This  area  is  not  empty  in  the 
actual  photography,  but  is  obscured  by  the  fireball  and  debris.  The  first 
time  any  fragment  is  observed  is  61  ms  after  time  zero.  This  supports  the 
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Figure  20.  Concrete  radial,  posttest  (Close-in). 


Figure  21.  Concrete  radial,  posttest  (Far-out). 
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Figure  23.  Posttest  view  of  upper  steel  column 
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TABLE  3.  GROUND-BASED  PHOTOGRAPHIC  DATA 


FRAGMENT 

APPARENT  AREA  (M^) 

FLIGHT  PROFILE 

AVERAGE  VELOCITY 

1 

0.1295 

STAB 

560 

2 

0.1403 

STAB 

629 

3 

0.1040 

TUMB 

474 

4 

0.1365 

TUMB 

381 

5 

0.0350 

TUMB 

164 

6 

0.0242 

TUMB 

359 

7 

0.0826 

TUMB 

590 

8 

0.0784 

TUMB 

382 

9 

0.0918 

STAB 

448 

in 

0.0252 

STAB 

371 

11 

0.1062 

ROTA 

423 

12 

0.1566 

STAB 

337 

13 

0.1113 

TUMB 

448 

14 

0.0085 

TUMB 

346 

15 

0.1056 

STAB 

276 

16 

0.1281 

ROTA 

402 

17 

0.1804 

STAB 

307 

18 

0.0656 

STAB 

279 

19 

0.0518 

TUMB 

281 

20 

0.1344 

TUMB 

172 

21 

0.0900 

88 

22 

0.0942 

TUMB 

138 

23 

0.0848 

TUMB 

62 

24 

0.1176 

TUMB 

152 

25 

•  0.0812 

STAB 

74 

26 

0.1825 

TUMB 

144 

27 

0.1020 

ROTA 

99 

28 

0.0384 

TUMB 

107 

29 

0.1030 

TUMB 

173 

30 

0.1692 

TUMB 

194 

31 

0.0768 

TUMB 

108 

32 

0.0468 

TUMB 

115 

33 

0.0630 

ROTA 

125 

34 

0.0330 

ROTA 

117 

35 

0.1219 

ROTA 

65 

36 

0.1178 

ROTA 

no 

37 

0.1320 

TUMB 

74 

38 

0.0344 

TUMB 

87 

39 

0.0504 

TUMB 

95 

40 

0.1197 

TUMB 

73 

41 

0.1176 

STAB 

71 

42 

0.0858 

ROTA 

51 

43 

0.1512 

TUMB 

510 

44 

0.0360 

TUMB 

359 

45 

0.0288 

TUMB 

347 

46 

0.0784 

TUMB 

279 

47 

0.1311 

TUMB 

393 

48 

0.1008 

TUMB 

379 

49 

0.1460 

TUMB 

322 

50 

0.0988 

354 
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TABLE  3. 

CONCLUDED 

FRAGMENT 

APPARENT  AR 

51 

0.1080 

52 

0.1380 

53 

0.0969 

54 

0.2574 

55 

0.1885 

56 

0.0644 

57 

0.0826 

58 

0.0812 

59 

0.1296 

60 

0.0954 

61 

0.2668 

62 

0.1290 

63 

0.2580 

64 

0.2208 

65 

0.1311 

66 

0.2296 

67 

0.4797 

68 

0.4949 

69 

0.2712 

70 

0.0841 

71 

0.0840 

72 

0.0896 

73 

0.1120 

74 

0.1914 

75 

0.0858 

76 

0.1848 

77 

0.1521 

78 

0.1890 

79 

0.0378 

80 

0.1024 

81 

0.1148 

82 

0.3040 

83 

0.1449 

84 

0.1104 

85 

0.0594 

86 

0.3562 

87 

0.0828 

88 

0.0480 

89 

0.0574 

90 

0.2511 

91 

0.2635 

92 

0.2988 

(M^)  FLIGHT  PROFILE 


TUMB 

TUMB 

TUMB 

TUMB 

TUMB 

TUMB 

TUMB 

TUMB 

TUMB 

TUMB 


TUMB 

TUMB 

TUMB 

TUMB 

TUMB 

TUMB 

TUMB 

TUMB 

TUMB 


TUMB 

STAB 

ROTA 

TUMB 

TUMB 

TUMB 

TUMB 


TUMB 

TUMB 

TUMB 


AVERAGE  VELOCITY(m/s) 


426 

466 

308 

340 

308 

231 

288 

419 

199 

262 

260 

138 

116 

99 

70 

144 

80 

94 

629 

435 

330 

333 

327 

336 

297 

187 

165 

158 

269 

134 

116 

165 

66 

58 

48 

75 

38 

34 

54 

32 

33 
88 
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Figure  24.  Fragment  trajectories  from  camera  F3386. 


Figure  25.  Fragment  trajectories  from  camera  F3382. 
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Figure  26.  Fragment  trajectories  from  camera  F3365. 

belief  that  the  metal  does  not  begin  to  break  up  and  move  until  after  a  sub¬ 
stantial  loading  takes  place.  Additionally,  at  this  specific  time,  the  frag¬ 
ment  is  first  observed  at  a  distance  approximately  two  HOBs  away  from  the 
charge.  The  last  time  any  fragment  is  seen  moving  vertically  is  at  422  ms. 
Other  fragment  motion  is  recorded,  but  these  fragments  have  either  reached 
their  apogee  and  are  returning  downward  or  are  moving  more  horizontally  than 
vertically.  The  limit  to  any  tracking,  as  shown  in  the  sketches,  is  bounded 
by  the  field  of  view  of  camera  and  lens. 

Another  attempt  was  made  to  use  the  ground-based  photographic  data  for  a 
more  precise  analysis.  The  data  from  the  three  cameras  were  combined  in  an 
attempt  to  triangulate  fragment  trajectories.  This  method  was  successful,  but 
for  very  few  fragments.  These  data  are  accurate,  but  have  a  large  scatter 
(this  indicates  the  extremely  large  scatter  in  the  data  using  one  camera). 

The  data  from  the  triangulated  analysis  are  presented  in  Appendix  E. 

Thirteen  fragments  were  commonly  observed  by  two  or  more  cameras. 
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Table  4  presents  a  summary  of  the  triangl ul ated  data.  Of  the  25  fragments 
triangulated  the  minimum  fragment  size  was  0.0252  m^  (0.2713  ft^)  and  the 
maximum  size  was  0.4949  m^  (5.33  ft^)  and  a  mean  of  0.1271  m^  (1.37  ft^).  The 
minimum  and  maximum  velocities  were  53  and  1638  m/s  (174  and  5374  ft/s), 
respectively  and  a  mean  of  491.13  m/s  (1611.4  ft/s). 

The  ground-based  photography  did  provide  useful  information,  but  the 
errors  associated  with  the  data  are  large.  Because  of  the  masking  and  dimen¬ 
sions  of  the  fireball  and  debris,  fragments  could  not  be  tracked  from  their 
points  of  origin.  Single  camera  coverage  is  not  adequate  enough  to  define  the 
flight  trajectory  of  any  fragment.  Multiple  camera  coverage  provides  more 
accurate  flight  trajectories,  however,  it  is  difficult  to  correlate  one  spe¬ 
cific  fragment  from  one  camera  to  other  cameras.  Also,  when  using  multiple  or 
a  single  camera,  it  is  uncertain  if  that  fragment  trajectory  is  in  one  plane. 

In  an  attempt  to  validate  the  theoretical  predictions,  the  photographic 
data  present  some  values  provided  the  assumptions  are  valid.  Using  the  time 
and  position  data  of  the  center  top  object  from  the  triangulated  data  a  decel¬ 
eration  rate  can  be  deduced.  For  simplicity,  and  lack  of  substantiating  data, 
assume  that  the  deceleration  is  linear;  also,  assuming  that  the  movement  of 
the  fragment  is  primarily  in  the  Z-direction,  and  horizontal  motion  is  negli¬ 
gible.  Extrapolating  the  measured  data  back  to  a  time  of  zero  and  subtracting 
the  calculated  deceleration  from  the  measured  velocity  it  is  found  that  the 
estimated  initial  velocity  is  878  m/s.  This  value  is  near  the  extrapolated 
information  from  the  Gurney  data.  It  is  important  to  realize  that  any 
fragment  could  be  used  to  estimate  the  initial  velocity,  and  different  values 
would  be  deduced.  However,  this  fragment  was  chosen  because  of  its  nearly 
vertical  trajectory  and  likelihood  of  high  quality  data  from  the  triangulation 
analysis. 

A  further  analysis  could  entail  extrapolation  of  the  fragment  back  to  its 
original  location  and  outward  to  impact  points.  Additional  correlation  could 
also  be  attempted  connecting  this  data  to  the  fragments  found  during  the 
ground  search.  At  present  the  data  from  the  aerial  photography  are  not  avail¬ 
able  for  correlation  and  analysis. 
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TABLE  4.  TRANGULATED  PHOTOGRAPHIC  DATA 


FRAGMENT 

APPARENT  AREA  (M^) 

AVERAGE  VELOCITY  (m/s) 

1 

0.1512 

572 

2 

0.0360 

402 

3 

0.0288 

369 

4 

0.1311 

420 

5 

0.1008 

315 

6 

0.0954 

361 

7 

0.2712 

8 

0.1104 

53 

9 

0.2296 

152 

10 

0.0826 

314 

11 

a 

1638 

12 

0.4949 

572 

13 

0.1062 

402 

14 

0.0252 

369 

15 

0.0784 

420 

16 

0.1062 

315 

17 

0.1113 

361 

18 

0.1403 

19 

0.0344 

53 

20 

0.1499 

152 

21 

0.1120 

314 

22 

a 

1638 

23 

0.1344 

152 

24 

0.0656 

314 

25 

a 

1638 

a  Undefined  area 
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Ground  Search  Ana1ysis--Inmediate1y  after  the  Pre-DIRECT  COURSE  event  an 
extensive  ground  search  was  initiated.  The  search  started  at  ground  zero  and 
extended  out  to  609  m,  and  consisted  of  several  men  walking  various  radials, 
tagging  fragments.  It  was  fairly  extensive  and  consisted  of  approximately  60 
manhours.  After  the  physical  ground  search  failed  to  yield  many  fragments  per 
manhour,  the  survey  process  began.  Measurements  of  angle  and  range  from 
ground  zero  were  recorded  within  an  accuracy  of  0.305  m  (1  ft)  in  range  and 
1  deg  of  angle. 

This  search  resulted  in  161  fragment  recoveries.  The  data  included  the 
remaining  portion  of  the  upper  steel  column  to  nuts  from  guy  cable  attachments. 
The  specific  number  of  types  of  fragments  is  listed  below; 

Number  Description 


31 

Angl  e 

47 

Channel 

51 

Metal  Plate  (could  be 

undefinable  channel 

or  angle) 

11 

Cable 

13 

Miscellaneous  (rings. 

nuts,  etc) 

8  Definite  portion  of  upper  steel  column 

The  fragments  ranged  from  ground  zero  out  to  518  m  (1700  ft).  A  detailed 
listing  of  the  fragments  recovered  is  presented  in  Appendix  F.  An  important 
consideration  is  that  all  fragments  were  not  recovered.  Many  were  believed  to 
be  buried  in  the  ground,  especially  on  the  dusty  radial.  Also,  approximately 
only  85  percent  of  the  total  mass  of  the  upper  steel  column  was  recovered. 

The  missing  15  percent  of  the  column  could  have  been  fragmented  and  recovered 
but  not  recognized,  or  simply  never  found. 

Figure  27  shows  a  plan  view  of  the  fragment  locations.  A  close  examina¬ 
tion  shows  that  a  preferential  scattering  took  place.  Four  areas  on  the 
figure  show  a  lower  distribution  than  the  remaining  areas.  These  areas  are 
45  deg,  135  deg,  225  deg  and  315  deg.  Review  of  the  initial  test  conditions 
reveals  that  these  angles  roughly  correspond  to  the  location  of  the  legs  of 
the  lower  battened  column. 
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Figure  27.  Plan  view  of  fragment  impacts. 
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The  range  distribution  of  fragments  was  also  calculated  via  a  Gaussian 
density.  The  results  are  shown  in  Figure  28.  The  mean  distance  of  fragment 
location  was  68.4  m  (224  ft)  from  ground  zero  with  a  standard  deviation  of 
11.3  m  (37.1  ft). 

Probability  of  Impact--Probabi 1 ities  of  impact  must  be  decided  to  determine 
if  hazards  to  personnel  or  experiments  exist  from  fragments.  As  stated  pre¬ 
viously,  161  fragments  were  recovered  from  the.  ground  search.  Assuming  that 
the  search  only  revealed  one-half  of  all  the  fragments,  the  total  count  would 
be  322.  By  dividing  the  test  bed  area  into  circles  one  can  determine  how  many 
fragments  landed  in  each  circle  (using  the  factor  of  safety  of  2,  a  conser¬ 
vative  estimate  should  be  formulated).  From  the  previous  section  a  scaled 
Gaussian  density  plot  provides  a  probability  of  occurrence  at  any  one  location 
based  on  all  of  the  data. 

Table  5  presents  the  area  data  and  probabilities  of  a  fragment  impacting 
in  that  area.  The  first  column  is  the  radius  of  the  area  of  concern.  The 
second  column  presents  the  normalized  Gaussian  density  from  Figure  28.  The 
next  column  presents  the  areas  from  the  radii  of  column  1  (A  =  iT(rQ^  -  r.j^)). 

The  fourth  column  shows  the  probability  of  a  fragment  impacting  in  the  area 
presented  in  the  third  column.  The  final  column  presents  the  probability  of  a 
fragment  impacting  in  the  corresponding  areas  assuming  the  ground  search  was 
only  one-half  complete.  According  to  this  type  of  analysis  the  highest  proba¬ 
bility  of  impact  is  between  the  60.9  and  121.9  m  (200  and  400  ft)  ground  range. 
Note  that  this  analysis  considers  the  number  of  fragments  which  land  in  the 
areas  and  not  the  probability  of  one  fragment  landing  in  the  area. 

If  the  potential  for  a  fragment  impacting  upon  a  certain  object  is  con¬ 
sidered,  a  further  step  is  required.  Knowing  the  object  area,  this  can  be 
used  as  the  numerator,  with  the  probability  from  the  right  column  of  Table  5 
as  the  denominator,  to  calculate  the  probability  of  impact  upon  the  object.  As 
an  example  consider  the  human  fragment  impact  hazard.  As  a  general  guideline 
the  safety  criteria  for  an  individual  susceptible  to  a  lethal  fragment  is 
defined  as  a  material  possessing  an  impact  energy  of  78.6  J  or  greater. 

An  individual  is  considered  to  occupy  0.56  m^  (6.03  ft^ ) .  Therefore,  not  more 
than  one  lethal  fragment  may  be  delivered  per  5.57  x  10  m  (6.0  x  10  f t  ) . 
Table  6  presents  the  probability  of  a  0.56  m  area  being  struck  by  a  lethal 
fragment . 

*0ppedahl ,  P.  E.,  Ltr  to  White  Sands  Missile  Range  (STEWS-NR-P) ,  N  Mex,  Sep  1982. 
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Figure  28.  Gaussian  density  fragment  distribution. 
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TABLE  6.  PROBABILITY  OF  HUMAN  HAZARDS 

range  probability  of  an  individual  in 

fml  PARTICULAR  RANGE  BEING  STRUCK 


0  - 

7.6 

1.9395  X 

10-^ 

7.6  - 

15.2 

1.2689  X 

10"® 

15.2  - 

22.8 

1.6293  X 

10-^ 

22.8  - 

30.5 

8.9634  X 

10-^  • 

30.5  - 

60.9 

5.2380  X 

10"^ 

60.9  - 

121.9 

1.8690  X 

10"^ 

121.9  - 

182.9 

1.5750  X 

10'“ 

182.9  - 

243.8 

8.9590  X 

10-^^ 

243.8  - 

304.8 

2.3377  X 

10"^^ 

304.8  - 

1828.7 

4.4806  X 

10-^^ 

Scaling  For  Direct  Course--The  simple  scaling  from  Pre-DIRECT  COURSE 
(2.177  X  10^  kg)  to  DIRECT  COURSE,  (5.445  x  10^  kg)  can  be  derived  by  Hopkinson 
or  cube  root  scaling.  This  implies  that  distances,  times,  and  impulses  are 
scaled  by  the  cube  root  of  the  ratio  of  the  energies  (yield).  Thus  the  cube 
root  scaling  factor  going  from  Pre-DIRECT  COURSE  scale  to  DIRECT  COURSE  scale 
is  2.924. 


The  scaling  factor  derived  above  is  valid  for  scaling  any  hydrodynamic 
variables  (p,  t,  and  I).  However,  it  is  not  appropriate  for  scaling  stresses 
and  strains.  These  quantities  follow  more  of  a  square  root  scaling  rule. 
Therefore,  stresses  in  the  steel  members  will  be  higher  for  DIRECT  COURSE 
(five  times). 

Even  though  the  stresses  are  higher  in  DIRECT  COURSE,  the  dimensions  of 
the  members  are  scaled  by  the  cube  root.  Therefore,  the  masses  are  also 
scaled.  The  major  differences  in  the  two  structures  will  be  along  welded 
joints.  At  these  locations,  depending  upon  the  design  criteria,  the  welds 
could  be  actually  stronger  than  the  steel  members.  Thus,  various  forms  of 
fragmentation  may  take  place. 
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If  we  look  at  the  basic  displacement  equation 
1  2 

s  =  “  at  +  V„t 
2  ° 

where  s  =  maximum  distance  traveled,  a  =  acceleration,  and  t  =  time.  Applying 
the  boundary  condition  that  at  t  =  0,  Vq  =  0  we  obtain 


2 


However,  the  acceleration  of  a  fragment  is  universally  proportional  to 
the  fragment  thickness.*  The  fragment  thickness  is  proportional  to  the  cube 
root  of  the  charge  weights,  as  is  the  time.  Therefore,  the  terms  cancel  with 
each  other  and  we  obtain  a  scale  factor  of  1,  which  indicates  that  the  veloc¬ 
ities  at  the  end  of  the  acceleration  phase  should  be  identical  for  Pre-DIRECT 
COURSE  and  DIRECT  COURSE. 

However,  since  the  fragment  thickness  is  proportional  to  the  cube  root  of 
the  charge  weight  we  again  derive  the  scale  factor  of  2.924.  (Since  the 
dimension  of  steel  thickness  were  scaled  from  Pre-DIRECT  COURSE  to  DIRECT 
COURSE  by  a  factor  of  2.924.)  Therefore,  the  ultimate  scaling  factor  for 
deriving  fragment  distances  for  DIRECT  COURSE  should  use  the  Pre-DIRECT  COURSE 
data  multiplied  by  2.924. 


*Ganong,  G.  E.,  (FCDNA/FCTOH) ,  Personal  Notes,  31  f4ay  83 
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IV.  CONCLUSIONS 

Fragment  prediction  and  analysis  is  unpredictable  and  almost  impossible 
given  one  test.  The  best  source  of  data  is  repetitive  testing  with  large 
amounts  of  information.  However,  using  various  technologies  and  sources  of 
information,  data  may  be  obtained  which  are  sufficient  for  rules  of  thumb 
analysis.  Quality  photographic  data  are  useful  but  may  be  too  advanced  con¬ 
sidering  the  state  of  the  art  of  fragment  predictions.  Simple  trajectory 
analysis  does  not  consider  the  drag  force  which  plays  a  major  role  in  fragment 
trajectories.  More  in-depth  analysis,  such  as  Reference  6,  provides  more 
accurate  results  even  when  extrapolated.  No  data  exist,  at  close  ranges, 
which  will  substantiate  the  assumed  initial  conditions.  If  such  data  did 
exist  a  more  valid  assessment  could  be  assumed. 

Notwithstanding  these  limitations,  enough  data  were  obtained  from  Pre- 
DIRECT  COURSE  to  give  a  factor  of  confidence  for  DIRECT  COURSE  as  currently 
envisioned.  Simply  stated,  if  a  person  is  not  within  a  121.9  m  radius  from 
ground  zero,  that  person  would  have  a  probability  of  less  than  1  in  63  billion 
of  being  hit  by  a  lethal  fragment. 
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APPENDIX  A 

This  appendix  presents  a  simplified  approach  to  the  trajectory  of  a  fragment. 
The  assumptions  used  are  important.  Although  this  report  is  not  based  on  this 
particular  analysis,  it  is  presented  because  of  its  common  occurrence. 
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SIMPLE  TRAJECTORY  ANALYSIS 

For  a  simplified  approach  to  trajectory  analysis  turn  to  the  equations  of 
motion  based  on  a  formulation  of  Newton's  second  law  of  motion. 

2  "^EXTERNAL  " 

since 

—  dv  — 

a  =  dt  = 

—  dr  .1 

v  =  ~  =  ^ 

a  =  r 

so 

^  "^EXTERNAL  =  "i*" 

Now  determine  the  forces  which  act  on  the  fragment.  The  net  force  will  be 
a  combination  of  aerodynamic  forces,  gravitational  attractions,  and  miscella¬ 
neous  effects  due  to  electromagnetic  fields,  thermal  radiation,  etc. 

Mathematically; 

^  ^EXTERNAL  "  "  ^1  i f t  ^drag  ^thrust  ^gravity  ^other 

Now  applying  the  following  assumptions; 

Fother  effect  of  miscellaneous  forces,  is  so  small  compareo  to 

the  other  effects  that  it  can  be  ignored  without  changing  the 
accuracy  of  the  solution. 

^lift  =  0  Since  the  fragments  are  assumed  to  have  an  average  angle  of 
attack  of  zero. 

^gravity  noticeable  gravitational  effect  is  due  to  the  attrac¬ 

tion  of  the  earth. 
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(■"gravity  =  "^9  =  constant 


where,  g  =  9.807  m/s^) 


■"thrust  ■  fragment  has  no  power  and  cannot  accelerate. 

Fdrag  =  1/2  P  Cj  A  (Vj/Vj) 


where , 


p 


A 


V 


a 


f^a/Val 


atmospheric  density 
coefficient  of  drag 

fragment  cross-sectional  area  as  viewed  from  along  the  flight 
path 

fragment  speed  with  respect  to  the  air  mass 

a  unit  vector  along  the  instantaneous  flight  path.  Drag  acts 
opposite  to  this 


If  the  fragment  is  unstable  (tumbles  or  rolls)  A  and  Cq  are  variables,  and 
■"drag  c*^3nges  randomly. 

Wind  is  not  a  force.  It  is  an  effect  wKich  appears  in  the  drag  term 
(specifically,  Vg). 

At  this  point  the  equation  of  motion  is: 


This  is  a  second  order,  nonlinear,  vector  differential  equation.  This 
equation  cannot  be  solved  in  closed  form  and  needs  further  simplification. 

Assume  that  drag  is  a  small  quantity  compared  to  gravity  (low  and  A). 
This  leaves  one  term  in  the  vector  equation  of  motion. 

mr  =  mg  or  r  =  g 

This  equation  states  that  the  acceleration  on  the  fragment  is  due  only  to 
the  effect  of  gravity.  However,  it  does  not  tell  the  position  and  correspond¬ 
ing  velocity  at  points  along  the  trajectory.  The  equation  must  be  integrated 
to  produce  equations  giving  position  and  velocity  as  functions  of  time. 
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For  further  analysis,  it  is  convenient  to  introduce  a  coordinate  system. 
An  orthogonal  system  composed  of  three  mutually  perpendicular  unit  vectors  is 
useful  for  trajectory  problems.  Assume  that  the  x-axis  is  horizontal  and 
pointed  downrange,  its  y-axis  is  also  horizontal  but  pointed  cross  range,  and 
its  z-axis  is  vertical  (Fig.  Al). 


Z 


DOWNRANGE 


Figure  Al.  Simple  trajectory  coordinates. 


To  establish  fragment  time  and  position  split  the  last  equation  into  three 
scalar  equations,  corresponding  to  the  coordinate  system. 


X  =  0 


z  =  -g 


t 


t 


t 


X  =  C 


z 


=  C3  -  gt 
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Constants  of  integration,  Ci ,  C2  and  C3  must  be  evaluated.  This  is  done 
by  checking  the  valves  for  x,  y  and  z  at  the  initial  condition  t  =  0. 


At  t  =  0: 

Xq  =  V^.  initial  velocity  in  x-di recti  on 
y^  =  0  assume  initially  no  cross  range  motion 
z  =  V,  initial  velocity  in  z-direction 

u 

therefore , 


integrating  again, 


Similarly  evaluating  €4,  C5  and  €5  we  find  that 

X  =  V^.t  +  Xi  y  =  0  z  =-  I  gt  +  v^.t  +  Zi 

Where  x^-  and  z^-  give  the  known  starting  fragment  position. 

The  equations  for  x  and  z  give  position  components  as  functions  of  time 
since  fragment  movement. 

The  equations  for  x  and  z  give  velocity  components,  also  as  a  function  of 
time . 

X  =  V  z  =  -gt  +  V7. 

These  equations  do  not  give  any  intuitive  feel  for  the  geometry  of  the  tra¬ 
jectory.  The  geometry  can  be  deduced  only  from  a  form  of  z  =  a  function  of  x. 

Solving  the  x-equation  for  t: 
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Substituting  into  the  z-equation: 

■  ■  •  -■  (^) 


+  z,- 


this  is  a  form  similar  to 
2 

z  =  ax  +  bx  +  c 

which  describes  a  parabola.  Therefore,  the  trajectory  is  parabolic. 

Figure  A2  depicts  a  typical  two-dimensional  trajectory.  The  following 
defines  the  variables  shown  in  the  figure. 

V,-  =  initial  velocity 

0^  =  initial  flight  path  angle,  measured  up  from  the  horizontal 
Vx.  *  V,  cos 
Vzi  =  Vi  sin  <&i 

Xi  and  Zi  =  coordinates  of  starting  point 

Xf  and  Xj:  =  coordinates  of  point  downrange  at  same  height  as  initial  point 
(Zf  =  Zi) 


Rs  =  Xf  -  Xi 


symmetric  range 


Figure  A2.  Simple  flight  profile. 
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Substituting  in  the  z  =  f(x)  equation: 


= 


Zf  - 


+  Zi 


V^sin*^ 

V^cos-s^ 


Solving  for  as  a  function  of  and  : 

(One  solution  is  R^  =  0,  which  is  trivial.) 


or 


^  Rj  RjSin®^ 

i  ' 


1  _  Rs 


_  1 


—  g  —  =  sin4>iCOS$i  =  —  sinZ®! 
2  \/2  1^2  ’ 


A 

R<.  =  —  sin2'ti 
3 


In  this  form  it  is  seen  at  once  that,  for  a  given  velocity  V^- ,  the  range 
is  a  maximum  when  the  sin2'5.j  is  a  maximum. 

Since  the  sin  has  its  maximum  value  of  unity  for  an  angle  of  90  deg,  the 
angle  above  will  be  45  deg.  Furthermore,  the  range  for  any  angle  any 
number  of  degrees  greater  than  45  deg  will  be  equal  to  that  for  an  equal 
number  of  degrees  less  than  45  deg. 

Two  other  values,  often  of  concern  to  trajectory  analysis,  are  T,  time  of 
flight,  and  maximum  altitude.  The  remainder  of  this  discussion  will 

derive  expressions  for  these  values. 
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Since  the  time  required  to  reach  the  highest  point  is  equal  to  the  time 
required  to  fall  the  same  distance  (no  drag),  the  formula  for  free-fall  may  be 
used.  The  equation  for  an  object  falling  from  rest  is: 

Vz  =  gt 

by  solving  for  t  and  substituting  for  V,  find: 

V,  V.-sin*.- 
t  =  -^  =  — ! - 1 

g  g 

Since  t  is  the  time  to  rise,  or  the  time  to  fall,  the  total  time  of  flight 
will  be  2t.  Therefore  the  time  of  flight  T  is: 


T  = 


2V,sin«>; 


To  find  the  maximum  trajectory  height,  is  proportional  to  the  square 

of  the  time  and  is  given  by: 

^ax  ■  9^ 

from  before,  =  gt 

=  2gZn,ax 

solving  for  obtain 


7  A 

^ax  2g 


Substituting  for  V,  we  have 


^max 


V^-sin*^ 
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APPENDIX  B 

This  appendix  presents  the  calculational  steps  from  Bishop's  work  as  presented 
in  Reference  6.  It  is  meant  to  allow  for  simple  calculations  using  the  appro¬ 
priate  steps.  Knowing  the  basics  of  the  problem,  the  range  of  the  fragment 
can  be  estimated. 
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BISHOP'S  ANALYSIS  (Ref.  6) 

The  procedure  for  calculating  the  maximum  range  of  a  projectile  can  be 
developed  simply  from  the  characteristic  distance  equation 

eqn  B-1  =  (0.167)(S^)(Cj^'^)^  (Pp/p^) 


where  p,  =  Average  value  of  the  air  density  (g/cc)  along  the  trajectory 


Pp  =  Density  of  projective  (g/cc) 


(Cp."  )  =  Average  value  of  Cq"^  along  the  trajectory,  where  Cp  is  the 

drag  coefficient  at  any  point  along  the  trajectory 
(dimensionless) 


S  =  Effective  thickness  (in)  of  the  projectile,  averaged  along 
^  the  trajectory 

in  conjunction  with  the  rotation  correction  factor 


K 


-1 


eqn  B-2  =  (1.57)(1  +  Re)"! 

^f 


where  =  Characteristic  distance  for  the  tumbling  rotational  case 

=  Characteristic  distance  which  corresponds  to  a  nontumbl¬ 
ing  plate 

Re  =  Thick  plate  shape  factor  (dimensionless)  =  fa^e  area 

Taking  average  values  along  the  trajectory  (subscript  £)  the  corrected  charac¬ 
teristic  distance  is  determined  from  equation  B-2. 
eqn  B-3  (K^'^)^  =  ( 1 -57)  ( Kfl  )a  (I+Re)"! 

in  which  (K^”^)a  is  the  corrected  value  of  (Kf"l)a  due  to  (tumbling)  rota¬ 
tion,  and  (Kf"l)a  is  defined  by  an  equation  similar  to  equation  B-1. 


eqn  B-4  (Kf-l)a  =  (0.167)  (Sf)(CD-l)a  {pjp) 

p  a 
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Equation  B-4  is  identical  to  equation  B-1  except  Sa  is  replaced  by  Sf.  The 
quantity  Sf  is  defined  as  the  fragment  (flat-plate)  volume  divided  by  the  face 
area  Af.  Thus  Sf  is  the  same  as  the  quantity  Sa  with  A  replaced  by  Af. 
Therefore,  Sf  is  constant  and  is  practially  equal  to  the  fragment  thickness  S 
(measured  in  inches),  unless  the  fragment  is  bent  or  curled. 

When  the  fragment  consists  of  a  bent,  curled,  or  folded  plate,  the  pre¬ 

vious  definitions  of  S  must  be  extended  in  order  to  estimate  the  corresponding 
range  distances.  The  effective  fragment  thickness  Sf  has  been  defined  as 
eqn  B-5  Sf  =  ^f/^f 

where 

Sf  =  Effective  thickness  (inches) 

Vf  =  Fragment  volume  (in^) 

Af  =  Face  area  (in2) 

For  a  flat  plate,  the  face  area  Af  is  unambiguous,  but  the  "face  area"  Af  of  a 
bent,  curled,  or  folded  plate  is  as  yet  undefined,  and  has  no  meaning.  In 
order  to  extend  the  definition,  we  refer  to  the  following  equation: 

eqn  B-6  A  =  Af  cos  H  +  Ae  sin  H 

for  0°  1  H  i  90° 

which  gives  the  frontal  projected  area  A  in  terms  of  the  two  areas  I\q  and  Af. 

From  equation  B-6,  it  is  clear  that  the  "edge-area"  A0  is  the  minimum 
frontal  projected  area,  and  this  minimum  occurs  when  the  inclination  angle  H 
is  equal  to  90  degrees.  This  conclusion  applies  as  yet  only  to  a  flat  plate. 
By  analogy,  let  us  define  Ag  as  the  minimum  frontal  projected  area  of  a  bent, 
curled,  or  folded  plate  and  let  us  choose  the  angle  H  so  that  H  is  90  degrees 
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when  the  bent,  curled,  or  folded  plate  presents  its  minimum  frontal  projected 
area.  Then,  in  order  to  preserve  the  form  of  equation  B-6,  we  define  Af  as 
the  frontal  projected  area  of  the  bent,  curled,  or  folded  plate  when  the 
inclination  angle  H  is  zero.  Obviously  this  occurs  when  the  fragment  has 
turned  90  degrees  from  the  position  where  it  presented  its  minimum  frontal 
projected  area,  defined  as  Ae.  It  follows  from  this  extended  definition  that 
the  generalized  "face  area"  Af  is  always  approximately  equal  to  the  maximum 
frontal  projected  area  of  the  fragment. 

Equation  B-4  gives  the  value  of  K"^  averaged  along  the  trajectory,  which 
results  in  equation  B-1  with  subscript  and  also  with  the  variable  projected 
area  A  replaced  by  the  constant  face  area  Af.  Hence  the  subscript  f  appears 
in  eqution  B-4  which  refers  to  a  nontumbling  fragment. 

Substituting  the  expression  for  (Kfl)  from  equation  B-4  into  B-3  we  ob¬ 
tain  the  characteristic  distance  (Kp'^a  for  ^  rotating  (tumbling)  projectile 

eqn  B-7  K  =  (0.262) ( Sf)(CD-l)a (  P/  a)(l+RE)-l 
eqn  B-8  K  =  (K^-l)a 

by  definition. 


The  general  method  for  calculating  a  typical  maximum  horizontal  range 
distance  (xr)  proceeds  according  to  the  following  schedule: 


Step  No. 

1  K 

2  r 

3  B(„ 


Calculation  Steps 

characteristic  distance,  feet,  from  equation  B-7  above 

from  r  =  K"1  (Vg^/g) 

from  Figure  B-1,  using  r  from  Step  2 
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(S33ii9aa) 

wa'39Nva  nmwixvw  aoj  3n9Nv  ivixini 
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FIGURE  B-l.  INITIAL  ANGLE 
FOR  MAXIMUM  RANGE 


AFWL-TR-83-52 


Step  No. 

4  Q 

5  w 

6  F 

7  (Zr) 

8  (xr) 


Calculation  (Continued) 

from  0  =  cos 

from  w  =  KQ 

from  F=2r  sin 

from  Figure  B2 

maximum  horizontal  range,  feet,  from  (xr)  =  (0.5)  w  (Zr) 


In  order  to  carry  out  all  the  steps  it  is  necessary  to  know  the  magnitudes 
of  the  following  quanitites: 


Sf  ^Face^area°^ ”  Effective  fragment  thickness  (inches) 

(CD)a  Df'^g  coefficient,  averaged  along  the  trajectory  (dimensionless) 

Pp  Projectile  density  (g/cc),  numerically  equal  to  the  specific 

gravity  relative  to  water 

Air  density  (g/cc),  averaged  along  the  trajectory 

Rc  Dimensionless  shape  factor.  Edge  area 

Face  area 

Vb  Initial  speed  of  fragment  (ft/sec) 

g  Acceleration  due  to  gravity 
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FIGURE  B-2.  SOLUTION  OF  RANGE 
EQUATION.  WIDE  RANGE  OF  Zp 
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This  appendix  contains 
All  the  variables  requi 
The  key  variables  are: 
(in),  thickness  (in)  of 
(ft/s),  and  distance  to 


APPENDIX  C 


the  results  of  calculations  using  Bishop's  methodology, 
red  (and  steps)  indicated  in  Appendix  B  are  presented. 
L,  W,  T,  VEL,  and  XR  which  refer  to  length  (in),  width 
the  fragment,  initial  velocity  of  the  fragment 
impact  point  (ft),  respectively. 
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APPENDIX  D 


This  appendix  contains  the  analysis  of  the  ground-based  photographic  data  as 
perceived  by  the  Physical  Sciences  Laboratory.  The  first  item  to  note  is  the 
camera  number  listed  on  the  second  line.  The  third  line  indicates  the 
fragment  number,  camera  number,  and  the  plane  angle  (respective  to  the  test 
bed  coordinates)  perpendicular  to  the  view  of  the  camera.  Also  presented  in 
the  fourth  line  is  the  deduced  photometric  size  of  the  fragment  and  its  flight 
attitude.  The  fragment  numbers  also  correspond  to  the  numbered  trajectories 
in  Figures  24  through  26.  The  data  present  the  position-time  and  calculated 
velocities  of  the  fragment  throughout  its  trajectory. 
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PRE-DIRLCT  tnUKSt  TEST,  7  JC T  19R?  3365 

FRAGMENT  TRAJECTURY  FROM  CAMERA  F-3365  DATA 
FRAGMENT  1  CAMERA  F-jj65  AZIMUTH  Gf  U-AXIS  3lZ.l  DEG 

APPARENT  SIZE:  L  =0.72N,  =0.21M  TUMtiLING 
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PRH-OIRtCT  CCURSc  TEST,  1  CC T  1902 
FRAGMENT  TKAJLCTGRY  FRUM  CAMERA  F-3365  DATA 
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AFWL-TR-83-52 


APPENDIX  E 

This  appendix  contains  the  analysis  of  the  triangulated  ground-based  pho¬ 
tographic  data  as  perceived  by  the  Physical  Sciences  Laboratory.  The  first 
line  on  the  data  page,  on  the  far  right  hand  side,  indicates  the  fragment 
letter  and  page  number.  The  second  line  presents  the  numbers  of  the  cameras 
from  which  data  were  obtained  for  this  analysis.  The  first  section  of  data 
gives  the  position-time  data  for  the  fragment  and  the  velocity  of  the  fragment 
corresponding  to  each  direction  of  motion.  Also  presented  is  the  calculated 
total  velocity,  angle  above  the  horizontal,  and  angle  with  respect  to  the  test 
bed  orientation.  The  second  set  of  data  (for  the  same  fragment  trajectory) 
lists  the  deviations  of  the  measurements  from  each  location,  the  error  of  the 
total  velocity,  and  the  locations  of  the  cameras  which  were  used  in  the  com¬ 
putation  . 
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AFWL-TR-83-52 


APPENDIX  F 

This  appendix  contains  the  results  of  the  ground  search  and  survey.  The  first 
column  identifies  the  fragment  number,  the  second  column  gives  the  azimuth  of 
the  fragment  location  (i  lO),  the  third  column  presents  the  ground  range  of 
the  recovered  fragment  (i  1  ft),  and  finally,  the  last  column  presents  general 
comments  on  the  description  of  the  recovered  fragment. 
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AZIMUTH  GROUND  RANGE 
ITEM  DEGREES  (FEET)  COMMENTS 


1 

1 

47 

2 

1 

54 

3 

5 

26 

4 

9 

16 

5 

9 

200 

6 

14 

46 

7 

15 

54 

8 

18 

29 

9 

18 

34 

10 

18 

124 

11 

19 

15 

12 

20 

24 

13 

29 

54 

14 

37 

36 

15 

37 

347 

16 

49 

89 

17 

49 

178 

18 

52 

29 

19 

54 

24 

20 

58 

215 

21 

58 

404 

22 

59 

277 

23 

63 

45 

24 

64 

120 

25 

64 

262 

26 

66 

■  25 

27 

71 

207 

28 

72 

23 

29 

80 

335 

30 

82 

152 

31 

84 

97 

32 

85 

118 

33 

86 

59 

34 

93 

309 

35 

94 

27 

36 

108 

231 

37 

112 

238 

38 

112 

307 

39 

112 

340 

40 

112 

416 

41 

113 

340 

42 

114 

61 

43 

118 

196 

44 

120 

18 

45 

125 

57 

46 

130 

117 

47 

136 

141 

48 

138 

95 

49 

141 

13 

50 

141 

142 

51 

142 

74 

52 

143 

242 

4x1  1/2**  angle 

20**  long  channel  **J1** 

21**  long  channel 
7x4**  steel  plate 
24**  long  channel  **L4** 

2x4**  steel  plate 

3x4**  steel  plate 

2x3  1/2**  angle 

4**  long  channel 

2x4**  steel  plate 

7**  long  angle 

1x2**  steel  plate 

2x1**  steel  plate 

20**  long  channel,  white 

**0**  ring  from  top  of  tower 

6**  long  channel  **K1** 

3x3**  steel  plate 
3x4**  angle 
1x3**  steel  plate 
1/2**  cable,  50  ft  long 
2x6**  channel 
2  ground  wire  clamps 
38**  long  angle 
48**  long  angle 
1  ground  wire  clamp 
6**  long  channel 
1x3**  steel  plate 
16**  long  channel  **A4** 

1x2**  steel  plate 

6x9**  lower  flange  plate 

24**  long  channel 

1x3**  steel  plate 

4x2**  steel  plate 

4x2**  channel 

5**  long  angle 

2**  ball  of  steel 

2x3**  steel  plate 

1x2**  steel  plate 

8**  long  channel  **C4** 

13**  long  channel 
2x2**  steel  plate 
9**  long  angle 
15**  long  channel  **J4** 

24**  long  channel 
4"  long  angle 
4x4**  steel  plate 
2x4"  steel  plate 
2x2"  steel  plate 
2x2"  steel  plate 
2x4"  steel  plate 
22**  long  channel  **J3** 

3/4"  cable  with  clamp,  1  ft  long 
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AZIMUTH  GROUND  RANGE 
ITEM  DEGREES  (FEET) 


53 

144 

231 

54 

145 

8 

55 

149 

106 

56 

152 

152 

57 

152 

578 

58 

153 

155 

59 

153 

165 

60 

153 

658 

61 

156 

7 

62 

156 

299 

63 

159 

125 

64 

160 

488 

65 

162 

109 

66 

167 

143 

67 

168 

10 

68 

169 

292 

69 

171 

323 

70 

171 

406 

71 

173 

120 

72 

173 

130 

73 

179 

45 

74 

180 

115 

75 

180 

435 

76 

180 

517 

77 

182 

310 

78 

183 

409 

79 

184 

417 

80 

187 

43 

81 

188 

454 

82 

189 

666 

83 

190 

43 

84 

190 

236 

85 

191 

44 

86 

193 

263 

87 

193 

486 

88 

196 

560 

89 

196 

781 

90 

198 

567 

91 

201 

9 

92 

201 

482 

93 

210 

1350 

94 

212 

627 

95 

213 

11 

96 

213 

13 

97 

213 

68 

98 

213 

73 

99 

213 

1700 

100 

214 

306 

101 

215 

320 

102 

215 

359 

103 

217 

491 

104 

222 

56 

COMMENTS 

6**  long  channel 
2x4**  steel  plate 
2x2**  steel  plate 
2x4**  steel  plate 
24**  long  copper  cable 
10  ft  long  copper  cable 
15  ft  long  copper  cable 
18**long  copper  cable 
3x2**  steel  plate 
4**  long  channel  **13** 

3x5**  steel  plate 

2**  long  angle 

8**  long  angle 

10**  long  channel 

7**  long  channel 

2x3**  steel  plate 

22**  long  channel,  white 

7**  long  channel 

24**  long  upper  pipe 

2x4**  steel 

3**  long  channel  **D2** 

6**  long  angle 

9**  long  channel 

20**  long  channel  **C3** 

3**  long  angle 

6**  long  side  off  channel 

7**  long  channel  **13** 

2x3**  angle 

3x2**  steel  plate 

16**  long  channel  **H3’* 

2x3**  angle 

4**  long  channel  **C3** 

2x4**  angle 
guy  cable  holder 
4x2**  steel  plate 
3/4**  steel  washer 
1x2**  channel 
9**  long  channel 
1x4**  steel  plate 
1x3**  steel  plate 
6**  long  channel  **3** 

9x5**  upper  pipe 
4**  long  plate 
5**  angle 
1x1**  angle 
2**  long  channel 
3/4-  nut 

10x10-  upper  steel  pipe 
2x5-  steel  plate 

24-  long  tether  cable  with  clamp 
2x3-  steel  plate 
7-  long  angle 
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